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As an atomic Bose Einstein condensate (BEC) is coupled to a source of 
uncondensed atoms at the same temperature and to a sink (extraction towards 
an atom laser) the idealized description in terms of a Gross-Pitaevsky equation 
(GP) no longer holds. Under suitable physical assumptions we show that the 
dissipative BEC obeys a Complex Ginzburg Landau equation (CGL) and for 
some parameter range it undergoes a space time patterning. As a consequence, 
the density of BEC atoms within the trap displays non trivial space time 
correlations, which can be detected by monitoring the density profile of the 
outgoing atom laser. The patterning condition requires a negative scattering 
length, as e.g. in ^Li. In such a case we expect a many domain collapsed 
regime, rather than a single one as reported for a closed BEC. 

I. INTRODUCTION 

The BEC dynamics in an atomic trap is ruled by a GP which in fact is a nonhnear 
Schrodinger equation (NLS) describing a conservative motion. Experimental evidence of 
BEC in a trap confirmed qualitatively a dynamical picture based on a GP description. 
On the other hand, extraction of BEC-atoms toward an atom laser introduces a dissipa- 
tion which must be compensated for by a transfer from the uncondensed fraction of trapped 



1 



atoms. Those ones on their turn must be refilled by a pumping process which, in the actual 
laboratory set ups is a discontinuous process but that we here consider as a continuous 
refilling, even though no working scheme is available yet. In Sec.^ we describe the addition 
of dissipative interactions through coupled rate equations, as done by Kneer et al. p. In 
Sec.^ we provide the physical grounds for an additional space dependent (diffusive) process 
and introduce an adiabatic elimination procedure, whereby we arrive at a closed equation 



which in fact is a CGL. In Sec.|3 we rescale the CGL around threshold for both positive 
(®^Rb) and negative (^Li) scattering lengths, showing that in the first case the BEG is sta- 
ble against space time variations, whereas in the second case the system can easily cross 
the instability barrier (so-called Benjamin- Feir line ||9|-[TT|). In Sec.0 we present numerical 
results showing that in the unstable case, rather than collapsing into one singular spot as in 
the isolated BEG, the open system presents many uncorrelated domains (space-time chaos). 
In Sec.^ we compare the strength of the nonlinear dissipative term introduced by us with 
the 3-body recombination rate. 

II. THE DYNAMICS OF AN OPEN BEC 

We know that a BEG is modeled by the GP 

ih^ = -^V'<j^ + V,.t<P + g\<j^\'<p. (1) 

where = ^/pe^^ is the macroscopic wave function describing the probability amplitude 
of the condensate, Vext is the trap potential, shaped as a harmonic oscillator with frequency 
u, and g is the coupling constant for the nonlinear (density | p dependent) self interaction. 
g is proportional to the s-wave scattering length 

g = tts. (2) 

m 

We discuss specific experimental situations concerning ®^Rb atoms (m = 1.4 10~^^ kg, 
= 5.77nm) and ^Li atoms (m = 0.11510^^^ kg, = — 1.45nm) [^. For an anisotropic 
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trap the frequency is = [uj^^Jyi^z)^ The above equation is formally a conservative 

non linear Schrodinger equation (NLS). Thus, it is a straightforward task to attach to a 
BEC all those space-time features familiar of a NLS as e.g. solitary structures and vortices 
||l^ , [T3|| , which have been explored in the recent past for a NLS, mainly in connection with 
pulse propagation in optical fibers ||14|,|15|. 



On the other hand the idealized picture of a BEC in an isolated system is in contrast 
with two physical facts, namely, 

i) The BEC is made of that fraction of atoms which have collapsed into the ground state 
(n = 0) of the harmonic oscillator trap potential; these atoms interact via collisions with 
those ones which are distributed over the excited states {n > 0) of the trap. The uncondensed 
atomic density Uu evolves in time not only because of the coupling with the condensed phase 
described by (p, but also because trapping and cooling processes imply a feeding (pumping) 
at a local rate R{f), the space dependence accounts for the non uniformities of the pumping 
process as well as for losses due to escape from the trap, at a rate 7„. 

ii) In order to have an atom laser, a radio frequency (rf) field is applied to the trap. The rf 
changes the magnetic quantum number of the atoms' ground state, thus transforming the 
trapping potential into a repulsive one and letting atoms escape from the BEC at a rate 7c 
ii- 

Both i) and ii) have been modeled by Kneer et al by adding dissipative terms to 
Eq.(|I]) and coupling the resulting equation with a rate equation for n^. In a slightly different 
formulation, this amounts to the following equations 

ih— = -— + Vext(p + g\ (p\^ (p- -h-fc(p + -^^nu4> (3) 
ot 2m 2 2 

and 

riu = R{r) - 7„n„ - Ln^nc. (4) 

here V is the rate constant coupling the condensed field with the uncondensed density n„, 
and ric =1 P is the local density of the condensed phase. We have modified the model of 



Ref. [§] as follows. At variance with [Q, where Eqs.(^,^ were written for the overall atomic 
population iV„ over the whole trap volume V, that is, 

iV« = y riud^r, (5) 

which is then coupled to 

Nc= f nj^r = f I |2 (r)dV, (6) 



here we prefer to deal with a local coupling. In fact Eqs.(|^J^) as written above are more 
convenient, as they refer to a local interaction. The coupling rate of uncondensed to con- 
densed atoms, r (m^ is given by the global rate used in which we call T', dividing 
by the trap volume V 

r' = ^(^-^). (7) 

Furthermore our local feeding rate R{r) is related to the overall rate Ru of ^ by 

Ru = J R{r)dh. (8) 

III. CGL PICTURE OF THE OPEN BEC 

Eqs.(^,^ were the basis of the model reported in Ref. [§. We wish to improve that 
picture, based on the following considerations. The uncondensed phase, n^^, is fed by a 
pumping process R{r) which is in general non uniform, and is locally depleted by its coupling 
with the condensed phase. As a result, n„ has a sensible space dependence and hence it 



undergoes diffusion processes. Precisely, by the fluctuation-dissipation theorem [|T6| , p!7| , the 
diffusion in velocity is given by 

I^. = 7n^(mV^). (9) 

m 

The corresponding diffusion constant in real space will be 
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= — = [m s ). (10) 

For ®^Rb at T = 100 iiK , and for 7^ = 5 10^s~^ (of the order of the average trap 
frequency Co ~ 7rl32Hz) this yields 

Dr = 210-'^{mh-^). (11) 

Thus we must add the term DrV^riu to Eq.(^. Once the BEG has been formed, the 
escape rate 7c in Eq.(|^) is compensated for by the feeding rate T Uu- As we set the BEG 
close to threshold, because of critical slowing down, the dynamics will be much slower 



than the n„ dynamics, thus we can apply an adiabatic elimination procedure find a 
quasi stationary solution for n„ in terms of and replace it into Eq.(^) which then becomes 
a closed equation for (p. 

We specify the above procedure by the following steps. First, rewrite Eq.(|) including 
diffusion 

riu = R{r) - -luriu + DrV'^Uu - Tuu | |^ • (12) 

Next, we take its space Fourier transform. The linear terms are trivial, whereas the nonlinear 
term should provide a convolution integral. However, since the condensate has collapsed on 
a single mode k = 0, then its Fourier transform is JF| 0(r) p =| |q 6{k) and Eq.(|T^ 
transforms as 

rik = Rk- (7« + Drk'^)nk - Tuk | |o • (13) 

where nk = J-'n^, Rk = J^R{r). The adiabatic elimination procedure consists in taking the 
stationary solution of Eq.(^) and replacing into Eq.(|). The stationary solution of Eq.(|l3|) 
is 

Uk = n r- (14) 

For long wavelength perturbations and far from saturation, the two additional terms in the 
denominator are less than unity. Here we consider a cylindrical volume with = Ly = 5/im, 



Lz = lO/iin (V = 250 (/iin)^) containing a condensate of Nc = 5 10'* atoms at a temperature 
T = lOOnK 0. It follows that k ^ 0{1/L) ^ 0.5 10^m~\ and hence Drk'^/^fu < 1. We can 
then expand Eq.(|l^ as 

n,^^il-^k'--\<l>\l). (15) 

7m 7m 7« 

The inverse Fourier transform of Eq. ([15|) is an operator relation as 

n. = -(l + ^V2--|0r). (16) 

7m 7m 7m 

As we replace this expression into Eq.(^), the operator acts on its right upon the space 
function (p. By doing this, we derive at a closed equation for (p which reads as 

where the square brackets contain the right hand side of Eq.(0). We now write explicitly 
the GP terms thus arriving at the following CGL 

2 7m ^ ^li 2m 27^ ^ 

The dissipative terms of Eq.(|TB|) represent respectively: i) difference between gain and 
losses, which implies a threshold condition; ii) a real diffusion which implies a spread of 
any local perturbation; iii) a real saturation term which provides a density dependent gain 
saturation. 



IV. RESCALED CGL NEAR THRESHOLD 



We herewith list the numerical values as taken from the experiment [jT9[ or from Ref. 
1^. We refer to a trap volume V = 0.25 10~^^m^; with a loss rate 7^ = 500 If we take 
the value |p F' ^ then F = T'V ~ 2 lO^^^m'^s"^. Furthermore, a reasonable estimate 

for the BEG escape rate toward the atom laser is |]^ 7c ~ 50 Therefore the threshold 
condition (gain=losses) is fulfilled for 
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RT/iu = 7c (19) 

where 12 10^^ m~^s~^, corresponding to Ru = RV ~ 3000 s~^. Finally we have for ^^Rb 
g/h = {4TTha,)/m = 0.48 m-^'^m^s~\ 

We can now write the parametrized CGL equation (|l^) in the dimensionless form 

= e0 + (1 + tci)V^(f) - (1 + ic2) I r 0. (20) 

where we have used the dimensionless time 

r = let, (21) 

and the dimensionless space coordinates 



X y z 
T' T' T' 



(22) 



where 



lo = {^Y" ^310-^m (23) 

is the characteristic length associated to the CGL dissipative dynamics, and the dimension- 
less condensated wave-function 

I^P^^I^P. (24) 
Note that tilde has been dropped in Eq. fpOl) . It follows from Eq.(|^) that 

1 , j?r . , , 

are the significant parameters of Eq. (PP]) . They are pure numbers. The term —iVext/^4> = 
—ico (f) in Eq. ([l8|) can be eliminated by a rotation transformation — > 06"*^^°*. 



We notice that Eg. (]20|) , derived by sound physical assumptions, is far from being a 
purely conservative (GP) or purely dissipative (real Ginzburg Landau) equation, but it 
displays both characters. 

However the Benjamin- Feir instability condition |p 

ci (-C2) > 1 (28) 

is not met by ^^Rb (ci=0.52, C2=2.05) and its dissipative CGL is fully inside the stable 
region. Hence the addition of dissipative terms may add interesting transient effects but 
does not lead to substantial qualitative changes with respect to the GP equation. Quite 
different is the case of ^Li (ci=0.04, C2=-0.51). Indeed even though the values of Ci and C2 
just listed give a stable dynamics, the fact that the scattering length is negative may lead to 
an instability if the parameters of the open BEG are slightly changed, e.g. if 7c is reduced 
by a factor 10 (which physically corresponds to a Li atom-laser with weaker losses). In such 
a case, we get Ci=0.4, C2=-5.1 and the open BEG is in the unstable region; we will denote 
this experimental situation by referring to an open Li* BEG. 



V. NUMERICAL SIMULATIONS 



As we have shown in Sec.|3, the coefficients of the GGL depend on the nature of the 
atoms forming the open BEG and also depend on the characteristic working parameters 
of the open BEG. Let us discuss the space time dynamics of the density of the condensed 
phase To do this, we proceed to the numerical integration of Eq . (|20D . The integration is 
performed on a two dimensional domain. This corresponds to a cross section of the 3-D cigar 
shape where the condensation takes place. This is justified by the fact that ujz» oJx ,ojy. 
The simulations are done on a 200x200 (Rb) or 256x256 (Li*) grid starting with an initial 
Gaussian distribution at the center of the domain. The numerical integration code is based 
on a semi-implicit scheme in time with finite difference in space. The chosen boundary 
conditions (at x = and y = ±Ly/2) are 
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dn 



0.10 



(29) 




where n is the normal at the boundary. Eq.(^) expresses the condition of an isotropic 

output flux of the condensed BEC (in the ideal situation of zero-gravity). The numerical 
coefficient on the right hand side of Eq.(^) is the dimensionless ratio between 7c and the 
velocity modulus of the condensed atoms, easily evaluated from the ground state solutions 



times for the Rb case (ci = 0.52, C2 = 2.05 and e = 0.5). The initial distribution evolves 
towards a stable quasi-uniform state. Fig.2 displays three cross section of Fig.l at different 
times, the solid line corresponds to the final stationary state and we observe the nearly 
uniform condensate on the overall domain. Fig. 3 illustrates a quite different situation: The 
values are now Ci = 0.4, C2 = —5.1 and e = 0.5 which corresponds to a Li* open BEC in 
the unstable region of use. The space-time chaotic dynamics emerges after a short transient 
(t < 10). Fig. 4 confirms that |0p is no longer symmetric with respect to x = (the same 
holds for the y-axis) . Figure 5 is aimed to show the spatial decorrelation of the signal when 
the condensate has entered the chaotic regime. 

The ID spatial power spectrum of the |0(a;)p function is shown in the lower curve of 
Fig. 5., the upper curve (solid line) is calculated by averaging the power spectrum of the 
function |0(x)p over a time interval from t = 500 until t = 1000 (taking a sampling time 
St = 10) within which the dynamics is statistically stationary. The results clearly indicates 
the large spatial decorrelation of the signal inside the chaotic regime. Indeed, it is well 

2 / 2 

known that the Fourier transform of a Gaussian function G oc '"^^ is again a Gaussian 
function G oc '""fc with afcCx = 1. On Fig. 5 it appears that the bandwidth in the Fourier 
space is much larger in the chaotic regime than for the initial distribution, which means a 
decorrelation of |0(x)p once the system becomes chaotic. 

To give a quantitative feeling, in the case of Li*, we have reduced 7c by a factor 10, which 
means that the normalization length is increased by vTO with respect to Eq. (pSf) , and it is 
/* ~ 10/im. In the numerical calculations we have considered a trap of linear size L = 50/*. 



of the harmonic oscillator 0. In Fig.l four snapshots of |0(x, ?/,t)p are shown at different 
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Since the ratio of the spectral widths between the chaotic and the initial spectra is about 4 
(estimated from Fig.5), it results that the coherence length in space time chaos is ~ 1/4 
the length of the initial Gaussian packet ctq- As seen from Fig. 4 ctq = 5/* = 1/lOL, hence 
(Tc — 1-2/* = 1/40L. These numerical estimates agree with a densitometric analysis of Fig.3. 



VI. DISCUSSION AND CONCLUSIONS 



Kagan et al. have discussed the collapse of a BEG in ^Li for a number of condensed 



atoms larger than the critical value Ncnt^ that is 



crit ~ 1 T- (30) 



This relation for Ncrit is obtained by equating the level spacing hu in the trap to the 
interparticle interaction energy no\g\ = N/V \g\ where g = Anti^as/m and V ~ af^^. Ref. ||20 



stabilizes the GP via a dissipative term corresponding to 3-body recombination processes. 
This amounts to a correction corresponding to a 6th power term in in a free energy 
potential. The dissipative equation of Ref. EO] is then 



(j) = —i[GP] — C,\(t)\^(t) + (pumping from the uncondensed portion). (31) 

In Eq. (|20|) , we have already treated the last term, here expressed in words, by the Kneer et 
al. approach Let us now compare the 5th power real damping entering Eq.(PTD with the 
3th power real damping of Eq. (pO|) . The cubic term is of the form G^cp, where 

G= = fSl^r = (32) 

Using the numerical values corresponding to the ^''Rb we obtain ~ 48 s^^ , for ^Li we 
have G3 ~ 1920 s^^. The cubic rate (G3) is a combination of the three characteristic rates of 
an open BEG. In a similar way we can introduce the rate G^ = C,\4'\'^- Taking the numerical 



values provided in Ref. |2^ we have the following ratios between the two dissipation rates 

G. I 480 for ^^Rb 



2 10^ for ^Li 
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(33) 



This result clearly indicates that for an open BEC, the 3-body recombination is negligi- 
ble with respect to the saturation cubic term that comes from the coupling between the 
condensed and uncondensed phase of the open BEC. 

To summarize, in this paper we have shown that in the framework of an atom-laser 
approach via two coupled equations, one for the uncondensed phase and the other one for 
the condensed phase, addition of a diffusion term for the uncondensed atoms and application 
of a proper adiabatic elimination procedure leads to a CGL dynamical equation for an open 
BEC. In the case of negative scattering length, a suitable adjustment of the escape rate 
implies entering the unstable regime of the CGL dynamics. Furthermore, within the chosen 
ranges of the parameters {e, 01,02), the 3-body recombination processes have a negligible 
influence. 

ACKNOWLEDGMENTS 

The Authors acknowledge E. Arimondo, T. Hansch, L. Pitaevskii and the BEC group 
at LENS Firenze for fruitful discussions. J.B. was partially supported by a EU Network 
grant (FMRXCT960010) "Nonlinear dynamics and statistical physics of spatially extended 
systems" and by the Belgian Programme on Interuniversity Poles of Attraction (PAI 04-6) 
initiated by the Belgian State Federal Office of Scientific, Technical and Cultural Affairs. 
This work was partially supported by INFM through the Advanced Research Project CAT. 



[1] E.P. Gross, Nuovo Cimento 20, 454 (1961); L.P. Pitaevskii, Zh.Eksp.Teor.Fiz. 40, 646 (1961) 
[Sov.Phys.JETP 13, 451 (1961)]. 

[2] F. Dalfovo, S. Giorgini, L.P. Pitaevskii and S. Stringari, Rev. Mod. Phys. 71, 463 (1999). 

[3] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman and E.A. Cornell, Science 269, 
198 (1995). 

11 



[4] C.C. Bradley, C.A. Sackett, J.J. Tollett and R.G. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 

[5] K.B. Davis, M.O. Mewes, M.R. Andrews, N.J. van Druten, D.S. Durfee, D.M. Kurn and W. 
Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

[6] M.O. Mewes, M.R. Andrews, D.M. Kurn, D.S. Durfee, C.G. Townsend and W.Ketterle, 
Phys.Rev.Lett. 78, 852 (1997). 

[7] I. Bloch, T.W. Hansch and T. Esslinger, Phys.Rev.Lett. 82, 3003 (1999). 

[8] B. Kneer, T. Wong, K. Vogel, W.P. Schleich and D.F. Walls, Phys.Rev. A58,4841 (1998). 

[9] T.B., Benjamin and J.E.Feir, J. Fluid Mech.27, 417 (1967). 

[10] B. Shraiman et ai, Physica (Amsterdam) D57, 241 (1992); H. Chate, Nonlinearity 7, 185 
(1994); H. Chate, in Spatiotemporal Patterns in Nonequilibrium Complex Systems, edited by 
P.E. Cladis and P. PalfFy-Muhoray (Addison- Wesley, New York, 1995). 

[11] H. Chate and P. Manneville, Physica A 224, 348 (1996). 

[12] R. Dum, J.I. Cirac, M. Lewenstein and P.Zoller, Phys.Rev.Lett. 80, 2972 (1998). 

[13] R.A. Sammut, A.V. Buryak and Yu.S. Kivshar, J.Opt.Soc.Am.B 15, 1488 (1998). 

[14] A. Hasegawaand and F. Tappert, Appl. Phys. Lett. 23, 142 (1973). 

[15] A. Kumar, Phys. Rep. 187, 63 (1990). 

[16] R.Kubo, J.Phys.Soc. Japan 12, 570 (1957). 

[17] T. Callen and T.A. Welton, Phys.Rev. 83 , 34 (1951). 

[18] H. Haken, Advanced Synergetics, Springer- Verlag, Berlin (1983). 

[19] BEC group at LENS, Firenze, private communication (1999). 

[20] Y. Kagan, A.E. Muryshev and G.V. Shlyapnikov, Phys. Rev. Lett. 81 (5), 993 (1998). 



12 



CAPTIONS 

Fig.l Open BEG of Rb: Time evolution of |0p using Eq.(^ with ci = 0.52, = 2.05 
and 6 = 0.5. The time step is (A^ = 0.01). The initial condition is a Gaussian. The figures 
are coded using a grey scale (white corresponds to the maximum value of 

Fig. 2 Guts of |0p at ?/ = 0, same parameters as for Fig.l. The dotted line corresponds 
to t=l, the dashed line is for t=10 and the solid line (t=500) corresponds to the final sta- 
tionary state. 

Fig. 3 Open BEG of Li* (escape parameters adjusted to be in the unstable regime): Time 
evolution of |0p using Eg. (pO|) with ci = 0.4, C2 = —5.1 and e = 0.5. 

Fig. 4 Guts of |0p at y = 0, same parameters as for Fig.3. The dotted line corresponds 
to t=l, the dashed line is for t=10 and the solid hne (t=500) corresponds to a chaotic state. 

Fig. 5 Spatial power spectrum of \4>{x,y = 0,t)p for Li*. The lower curve (dotted line) 
corresponds to t=l. The upper curve (solid line) is calculated by averaging the spatial power 
spectrum of \(p{x,y = 0,t)p between t=500 and t=1000 (by steps of 6t = 10). 



13 



figure i : r.L. Areccm, J. rJragara ana L. Lyasteiiano. 
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Figure 2 : F.T. Arecchi, J. Bragard and L. Castellano. 
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figure d : r.L. Areccm, J. rJragara ana L. Lyasteiiano. 
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Figure 4 : F.T. Arecchi, J. Bragard and L. Castellano. 




Figure 5 : F.T. Arecchi, J. Bragard and L. Castellano. 




